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Abstract The properties of superdense matter in neu- 
tron star (NS) cores control NS thermal states by affect- 
ing the efficiency of neutrino emission from NS interi- 
ors. To probe these properties we confront the theory of 
thermal evolution of NSs with observations of their ther- 
mal radiation. Our observational basis includes cooling 
isolated NSs (INSs) and NSs in quiescent states of soft 
X-ray transients (SXTs). We find that the data on SXTs 
support the conclusions obtained from the analysis of 
INSs: strong proton superfluidity with T^'^^ > 10^ K 
should be present, while mild neutron superfluidity with 
j^max ^ 2 X (lO^-lO^) K is ruled out in the outer NS 
core. Here T^'^^ and T^'^'^ are the maximum values of 
the density dependent critical temperatures of neutrons 
and protons. The data on SXTs suggest also that: (i) 
cooling of massive NSs is enhanced by neutrino emission 
more powerful than the emission due to Cooper pairing 
of neutrons; (ii) mild neutron superfluidity, if available, 
might be present only in inner cores of massive NSs. In 
the latter case SXTs would exhibit dichotomy, i.e. very 
similar SXTs may evolve to very different thermal states. 
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1 Introduction 

Neutron stars are very compact; their cores contain mat- 
ter with density p a few times larger than the standard 

This work is supported by Polish MEiN (grant no. 1P03D- 
008-27), by Russian Foundation for Basic Research (grants 
05-02-16245, 05-02-22003), by the Russian Federal Agency 
for Science and Innovations (grant NSh 9879.2006.2), and by 
Russian Science Support Foundation. 

K.P.Levenfish 

lofi^e Physical Technical Institute, St. -Petersburg, Russia 
E-mail: ksen@astro.ioffe.ru 

P.Haensel 

N.Copernicus Astronomical Center, Warsaw, Poland 
E-mail: haensel@camk.edu.pl 



nuclear matter density po = 2.8 x 10^** g cm~^. Many 
properties of this matter cannot be calculated precisely 
or studied in laboratory experiments. However, these 
properties can be constrained by co mparing neu tron star 
theory wit h observations; se e e.g., lYakovlev~fc Pcthic^ 
(|2004l ) and lPage etld] (|2006[ ). for recent reviews. 



In this paper, we describe current contraints on com- 
position and superfluidity of neutron star cores, which 
can be obtained by comparing calculated thermal states 
of neutron stars with observations of thermal radiation 
from middle-aged INSs and NSs in SXTs in quiescent 
states. INSs are thought to cool gradually from initial 
hot states via neutrino emission from NS interiors (at age 
< 10^ years) and via surface photon emission (at the el- 
der age; the so called neutrino and photon cooling stages, 
respectively). NSs in SXTs will be assumed to support 
their warm states owing to deep crustal heating (pyc- 
nonuclear reactions) in accreted matter fsee lBrown et all 
|1998[ ). Their thermal energy is partly emitted by neutri- 
nos and partly by the surface radiation. In both cases, 
thermal states of NSs are very sensitive to composition 
and superfluidty of matter in their cores. Although INSs 
and SXTs are different objects, their observations al- 
low one to test the same physics of superdense matter 
(Yakovlev et al. 2003; 2004). 

Composition and superfluidity of baryons in NS cores 
are main regulators of NS thermal states. These prop- 
erties are thought to be the same for a given density 
in all NS cores. The composition determines dominant 
processes of neutrino cooling in NSs. Superfluidity, if it 
appears, reduces rates of neutrino processes at work, but 
opens an additional neutrino mechanism associated with 
Cooper pairing of baryons. Superfluidity affects also the 
NS heat capacity. Other regulators of NSs thermal states 
(composition of NS heat-blanketing envelopes, strength 
and geometry of NS magnetic field, etc.) are also im- 
portant but, taken alone, do not allow one to recon- 
cile theory of NS thermal sta te with observations (see, 
e.g- lYakovlev k Pethicd2004f ) . Moreover, they may vary 
from one star to another. We will neglect them, for clar- 
ity. 
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Fig. 1 Thermal states of nonsuperfluid NSs compared to observations. Left: effective surface temperatures of INSs, redshifted 
for a distant observer, versus their age. Right: redshifted photon luminosities of NSs in SXTs in quiescence versus time- 
averaged mass accretion rate. The dotted curves refer to the basic NS model (a non-superfluid low-mass NS which cools 
slowly through the modified Urea process). Three solid curves on each panel display scenarios with the enhanced neutrino 
cooling (maximum-mass NSs with inner cores containing - from top to bottom - kaon condensates, pion condensates, and 
nucleons with large proton fraction, sufficient to open direct Urea process). Hatched regions between the basic curve and any 
solid curve can be filled by curves of NSs with different masses, from ~ 1 Mq to the maximum one, for a corresponding NS 
composition. 



2 Composition of NS cores 

An NS core can be divided into the outer core (p < a 
fewpo) a-iid the inner one (higher p). The outer core con- 
sists mostly of neutrons, with a small (a few %) admix- 
ture of protons and leptons. Low-mass NSs have only the 
outer core. In the absence of nucleon superfluidity, they 
cool down via the modified Urea processes and weaker 
neutrino processes of nucleon-nucleon bremsstrahlung. 

More massive NSs possess also the inner core whose 
composition is largely unknown. According to different 
hypotheses, the inner core may contain nucleons, hyper- 
ons, pion or kaon condensates, or quarks. Even if com- 
posed of nucleons and leptons, superdense matter can 
have a large fraction of protons. In each of these cases, 
neutrino emission processes of direct Urea type become 
allowed, much more powerful than the basic modified 
Urea processes. The most powerful is the nucleon (or 
hyperon) direct Urea process. Neutrino reactions with 
pions (kaons) are about two (four) orders of magnitude 
weaker. 

In the presented figures we compare various theo- 
retical predictions of thermal states of INSs and NSs 
in SXTs with observations. For INSs, we plot the red- 
shifted effective surface temperature versus NS age 
t. For SXTs, we show the redshifted surface luminos- 
ity Lg° of NSs in quiescence versus time-averaged mass 
accretion rate M. We assume that NSs in SXTs are 



in thermal steady-states, with the deep crustal heating 
(for a correspond ing M) balanced by n eutrino and pho- 
ton emission; see lYakovlev et al.l (|2003 ^ for details. The 
deep crustal heating is calculated using the model of 
iHaensel fc ZdunikI (Il99nh . In the majority of cases NS in- 
teriors are nearly isothermal owing to high thermal con- 
ductivity, with the main temperature gradients located 
in the heat-blanketing envelope near the NS surface. We 
omit technical details because of space restrictions. Ob- 
servat ions of INSs are the same as in iKaminker et al.l 
(I2002D . The SXTs Aql X-1, 4U 1608-52, SA X 1808 .3-3658 
and Cen X-4 are described in ,Yakovlev et al. I (12003.) . 
The data on XTE J2123-058, KS 1731-260, RX J1709- 
2639, SAX 1810.8-260 9 and IH 1905 -HOOO are taken from 
iTomsick et al. ( 20041) . ICackett et al.l (|I006i), Jonker et al. 
(2003, 2004 and 2006), respectively 



Fig.[T]shows thermal states of nonsuperfluid NSs with 
three different compositions in the inner core. In both 
panels, the upper (basic) dotted curve corresponds to 
low-mass NSs which possess no inner core and cool slowly 
via the modified Urea process from the outer core. The 
basic curve is almost independent of NS mass M as long 
as the inner core is absent. For higher M, we obtain 
noticeably colder NSs, cooling via enhanced neutrino 
emission from the inner core. Their surface temperature 
strongly depends on the composition in the inner core. 
The coldest is the maximum- mass NS (solid lines). 
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Fig. 2 Tiiermal states of nonsuperfluid INSs and SXTs with nucleon cores based on an EOS of Prakash et al. (1988). This 
EOS opens the nucleon direct Urea process at p > pD = 2.8po, i-e., in inner cores of NSs with M > Mo = 1.358 Mq. The 
thick dotted basic curves are the same as in Fig. [T] From top to bottom, thin solid curves correspond to NSs with masses 
growing from 1 Mq to Mmax = 1.977 M©, with a step of 0.02 Af0. An upper dense bundle of curves in each panel contains 
low-mass stars with M < Mu, while a less dense bottom bundle contains NSs with M > 1.06 AIu- 



The coldest INSs observed to date are consistent with 
all th ree neutrino emission scenarios (see lYakovlev et all 
l2003l for details of the models). At present, the data 
on SXTs contain colder sources and seem to be more 
restrictive. From the upper l imits on the thermal lumi- 
nosity of SAX 1808.4-3658 (ICampana et al.l [20021) and 
IH 1905-^000 (jJonker et al.ll2006t) . one can infer that a 
dominant process in superdense matter should be more 
powerful than direct-Urca-type processes with kaons or 
pions. However, these results should be taken with cau- 
tion - some issues of theory and observations of SXTs 
still have to be clarified (Sect. [5]). 

In what follows we will limit ourselves to the simplest 
nucleon models of NS cores descibed, e.g., in Yakovlev 
& Pethick (2004). Results presented in Figs. [2HH] are ob- 
tained using our generally relativistic code of NS thermal 
evolution. We assume (for clarity of our study) the ab- 
sence of light-element accreted envelopes on NS surfaces 
and neglect the effects of NS magnetic fields. NS models 
in Figs. BHU are based on a m oderately stiff EOS pro- 
posed bv lPrakash et al.l (|l988h . This EOS opens the nu- 
cleon direct Urea process at p > 2.8pQ, i.e., in the inner 
cores of massive NSs ith M > Md = 1.358 Mq. 

According to Fig. [21 models of non-superfluid NSs 
with enhanced emission in the inner cores cannot ex- 
plain the data on INSs and SXTs. First, they are un- 
able to interpret hottest sources. Second, a transition 
between widely spaced hot and cold NSs thermal states 
occurs within an unrealistically narrow NS mass range 
~ 0.01 A/q. As we show below, including superfluidity 
relieves these shortcomings. 



3 Superfluidity of nucleon matter 

Microscopic theories of dense nucleon matter predict that 
below some critical temperature, neutrons and protons 
in NS cores are superfluid. However, the critical temper- 
atures Ten and Tcj, as a function of density are very un- 



certain; see, e.g., iLombardo fc Schulzd (|200lh . Therefore, 
at present, it seems reasonable to rely on a few general 
points of recent theories of nucleon superfluidity: 

— proton pairing occurs in the state and persists 
from ~ 0.5po up to a few po! 

— critical temperature for protons, Tcp{p), may be rather 
high, with the maximum > 10^ K somewhere be- 
tween po and 2po; 

— neutrons form pairs in a ■^P2 state; this pairing is, 
typically, weaker and persists to higher densities than 
the proton one; 

— critical temperature for neutrons, Tcn{p), has maxi- 
mum somewhere between po and few po; as a rule, 
this maximum is shifted toward higher p relatively to 
the maximum of Tcp{p). 

These features of Tcp{p) and Tcn(p) can be simulated 
with phenomenological models. We adopt the models 
shown in Fig. [T] 

Strong p roton superfluidity. According to several au- 
thors (see lKaminker et al.ll2002l for references) a nucleon 
NS model with the open direct Urea in the inner core 
and strong proton pairing in the outer core can explain 
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Fig. 3 Same as in Fig. [2]but for strong proton superfluidity in the NS core (model "Ip" in Fig. (Tjl. In this case, Tcp{p) has 
maximum ~ 6.9 x 10^ K at p ~ 2 po, remains ^ 5 x 10* K at p Si 3.2 po, and dies out at p ~ 3.3 po- Accordingly, the entire 
cores of NSs with M ^ 1.49 are strongly superfluid, while the inner cores of NSs with M k, 1.52 have nonsuperfluid 
central kernels. 



available observations of INSs. Here we show that this 
model can explain also the data on SXTs. Our results 
are displayed in Fig. [31 

The effect of strong proton superfluidity is twofold. 
First, this superfluidity suppresses neutrino emission from 
NSs, making them hotter at a given age or mass accre- 
tion rate. This brings thermal states of slowly cooling 
low-mass NSs into the agreement with observations of 
hotter sources. In massive NSs with enhanced cooling, 
strong proton superfluidity may spread out the opening 
of the direct Urea process over some density range. As 
a result, the enhanced cooling sets in gradually with the 
growth of NS mass (not as sharp at M = Md, as in 
Fig. [2]) . This allows one to interpret the colder sources 
as massive NSs of different masses. 

Strong proton superfluidity, with Tcp > 7 x 10^ K, ap- 
pears in hot NSs with large neutrino luminosity, so that 
the neutrino emission due to Cooper pairing of protons is 
unimportant (Yakovlev et al. 2001). However, this sup- 
erfluidity suppresses the powerful Urea processes. Now, 
low-mass NSs mainly cool via a much weaker neutron- 
neutron bremsstrahlung (unaffected by proton pairing). 
It slows down the cooling and let the low-mass stars be 
hotter at a given <: or M. In this way one can inter- 
pret observations of the hotter sources without invoking 
any reheating mechanism. As noted by Yakovlev et al. 
(2004), the presence of light elements (H, He) in the NS 
heat-blanketing envelope facilitates interpretation of the 
hotter INSs. Such an envelope is more heat transparent 
and let the NS look hotter for a given inner temperature. 
Moreover, in that case even a weaker superfluidity (with 



Tcp <L 10^ K) allows us to interpret the hotter sources; 
we have checked that this is also true for SXTs. 

Thermal evolution of slowly cooling low-mass NSs is 
almost independent of the assumed EOS (Page & Ap- 
plegate 1992), as well of the NSs mass. At a given age or 
mass accretion rate, all low-mass stars have nearly the 
same inner temperature. Consequently, the appearance 
of superfluidity affects all these stars in the same way; 
cf. Figs. [Hand [3] (and Fig. [Jj below). This property holds 
as long as proton pairing is strong in the entire NS core. 

The situation is different in more massive NSs. The 
impact of proton superfluidity on these stars depends 
also on how far proton pairing extends into the inner 
core and how steep is t he slope of the TrJp) pr ofile 
in this superfluid region ( Yakovlev fc Haensel 20031 ). In 
Figs. [3| and [4| the inner core of NS models with Md < 
M < Mmax may occupy densities from 2.8po to 9.2po 
being superfluid at ^ (2.8-3.3)po and normal at higher 
p. Accordingly, in medium-mass NSs with M ^ (1.36- 
1.52) Mq the inner core is entirely superfluid. Superflu- 
idity is mild near the outer core boundary, and weakens 
rapidly toward the core center. The more massive the 
star, the weaker proton superfluidity in its central part 
and the less it suppresses the direct Urea process, thus 
letting the star cool down faster. In this way thermal 
states of medium-mas s INSs become dependent of M 
([Kaminker et al.ll200lh . 



In the innermost central kernels of most massive NSs 
(M > 1.52 Af©) proton superfluidity dies out, and the 
direct Urea process reaches its full power. This powerful 
fast cooling, even in a small part of the core, renders the 
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Fig. 4 Same as in Fig. |3]but with an addition of mild neutron superfluidity in the outer NS core (model "Int" in Fig. 
[7|. In this case, Tcn{p) has maximum at p ~2po and dies out at ~5.7po (i-e., neutrons are nonsuperfluid in the central 
kernels of NSs with M > 1.89 Af©). The maximum of the Tcn{p) curve is wide - Tcn{p) remains within (2.0-3.3) x 10* K 
at p ~ (0.72-3.2)po- This density range corresponds to central densities of low-mass NSs (M < 1.358 Mq) and moderately 
massive NSs (with M (1.36-1.48) Mq). 



superfluid NS as cold as the nonsuperfluid one; of. Figs. [5] 
and [21 



Mild neutron superfluidity in the outer NS core. 

Let us now add mild neutron superfluidity and assume 
first that this superfluidity is located in the outer core. 
We have adopted model "Int" shown in Fig. [71 Cooling 
of INSs with superfluiditi es "Int" and "Ip" was studied 
bv lKaminker et all (|200l[ ) and was shown to be inconsis- 
tent with hotter sources. We have tested this statement 
against the data on SXTs with the same conclusion. Re- 
sults are displayed in Fig. [H 

In cooling INSs the effect of mild neutron superfluid- 
ity is very spectacular. First, neutron pairing reduces the 
heat capacity of NSs, because the heat is stored mostly 
in neutrons. The reduction amounts to a factor of sev- 
eral in low-mass NSs. This makes low-mass NSs very 
cold after appearing of such pairing at t > 3 x 10* yr. 
The cooling is additionally accelerated by powerful neu- 
trino emission, triggered by Cooper pairing of neutrons. 
In result, the low-mass NSs models become unable to in- 
terpret four old sources: Geminga, RX J1856-3754, PSR 
1055-52 and RX J0720-3125. 

In the most massive NSs with M > 1.54 M© the effect 
of neutron superfluidity is opposite. In these stars the 
mild pairing appears earlier, at t ~ 10^ yr, when neutrino 
emission due to the direct Urea process is much stronger 
than the Cooper pairing neutrino emission (which is, 
therefore, insigniflcant). On the other hand, additional 
suppression of the direct Urea process by neutron super- 



fluidity makes the most massive NSs a bit hotter in the 
advanced stage of neutrino cooling era. 

As for transiently accreting NSs in SXTs, their ther- 
mal states become independent of the NS heat capacity 
as soon as they reach the steady-state regime (Yakovlev 
et al. 2003; 2004b). Hence, a strong reduction of the 
heat capacity does not affect NS thermal states at the 
photon-dominated stage. That is why the slope of the 
dotted curve at M<10~^^ Mq ■ yi^-^, corresponding to 
this stage, remains unchanged. However, at the neutrino 
stage, neutrino emission due to Cooper pairing of neu- 
trons makes low-mass SXTs much colder (which unables 
one to interpret such hot sources as Aql X-1 and RX 
1709-2639). The same effect in lightest among the mas- 
sive SXTs hampers the explanation of KS 1721-260. Sim- 
ilarly to high-mass INSs, high-mass SXTs become a bit 
hotter; neutron superfluidity, in its turn, slightly broad- 
ens the span of their allowed thermal states, stretched 
before by the proton pairing. 



Mild neutron superfluidity in the inner NS core. 

The presence of mild neutron superfluidity in the outer 
NS core contradicts observations of INSs and SXTs, as 
it renders low-mass NSs too cold. However, mild pairing 
might occur in the inner core which is present only in 
massive NSs. Let us explore this case. 

From Fig. [H we see that neutrino emission due to 
Cooper pairing of neutrons may initiate a very fast cool- 
ing and, thus, may serve as a fast cooling agent (at least, 
for INSs). 
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Fig. 5 Thermal states of superfluid NSs based on the EOS of lDouchin fc Haensell (120011 ) which prohibits the direct Urea 
process in NS cores. The figure illustrates the so called "minimal cooling scenario" in which Cooper pairing of neutrons 
operates as enhanced cooling agent. The scenario assumes strong proton pairing in the outer NS core and mild neutron 
pairing in the inner core (models "pi" and "ntl" in Fig. [T] respectively). In both panels the dotted basic curve refers to the 
same cooling scenario as in Fig. [T] From top to botton, thin solid curves show NSs models with masses growing from 1 Mq 
to Mmax ~ 2.051 Mq, with a step of 0.02 Mq. Pairing "pi" has a wide maximum at p ~ 1.6po and dies out at ~ 3.8 po, i.e. 
in the central kernels of NSs with M k, 1A9Mq. Strong pairing with Tcp ~ (2.0-6.9) x lO' K persists up to p ~3.2po and 
extends over entire cores of NSs with M ^ 1.27 M0. Neutron pairing of "ntl" is weak in the outer core, at p ^ 3.1 po (i.e., 
in NSs with M ^ I.24M0), has a sharp maximum of mild strength in the inner core, at p ~4.7po, and dies out rapidly 
as p approaches ~6.7po. Mild pairing with Ten ~ (2.0-6.0) x 10* K extends over the density interval ~ (3.1-6.4)po, which 
corresponds to the inner cores of NSs with M ~ (1.22-1.95)M0 or to spherical layers around the nonsuperfluid central kernels 
of NSs with M = (1.97-2.O5)M0. 



In INSs , such a possibil i ty wa s studied by I Page et al.l 
(|2004D and busakov et all (|2004l ). In particular, the lat- 
ter aut hors consider nucleqn NS models based on the 
EOS of lDouchin fc Haensell (|200lh which prohibits the 
direct Urea process. Their cooling scenario is based on 
strong proton pairing in the outer NS core and mild neu- 
tron pairing in the inner core (phenomenological mod- 
els "pi" and "ntl", respectively; see Fig. |7|). The model 
"ntl" has a specific dependence Tc„(p) which keeps neu- 
tron pairing weak in the outer core. Let us remind that 
weak superfluidity, with Ten ^ 2 x 10^ K, docs not ap- 
pear in low-mass NSs in the neutrino cooling era: these 
NSs are too hot. Accordingly, in this era neutron pairing 
"ntl" affects only massive NSs. 



iGusakoy et al.l ( 2004 ) compared such cooling scenario 
(called by Dany Page a "minimal cooling model" , for its 
simplicity, and the minimal number of its ingredients) 
with observations of INSs, and found it marginally con- 
sistent with them. The scenario can explain those INSs 
whose thermal emission is detected with confidence, as 
illustrated in the left panel of Fig.O However, these au- 
thors noted that the model will fail if the two sources, RX 
J0205-f 64 and CTA 1 (expected to be thermally emitting 
INSs), turn out to be much colder than their presently 



established upper limits, or if very cold INSs are detected 
in future. 

Since NSs in SXTs are expected to be rather mas- 
sive, owing to accretion of matter from their compan- 
ions, it is worth to test the "minimal cooling scenario" 
against the data on SXTs. That is done in the right 
panel of Fig. [51 One can see that the model success- 
fully reproduces hotter SXTs, including the frequently 
bursting sources AqlX-1, RX 1709-2639 and 4U 1608- 
52. However, the upper limits for at least two sources, 
IH 1905-hOOO and SAX 1808.4-3658, definitely fall far 
below the predictions of the model. 

The transient sou rce IH 1905-1-000 w as recently ob- 
served with Chandra ( Jonker et al.l200(t ): the upper hmit 
for its quiescent thermal luminosity seems to be firm. On 
the contrary, the value of M for this object is quite un- 
certain; some physical arguments allow one to believe it 
to be higher than 10~^^Mq • yr"-'^. For the frequently 
bursting transient SAX 1808.4-3658 the mass accretion 
rate is known more accurately, while the upper limit of 
its quiescent thermal luminos ity is less certain . It va ries 
from ^ 6 X 10^^ crg-s^^ in ICampana et al.l ( 20021 ) to 



~ 4 X 10 erg-s in estimations of P. Shtykovski (pri- 
vat e communication; see deta ils in lYakovlev et al" 1 120041 
and lYakovlev fc Pethic^l2004^ . Nevertheless, both these 
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Fig. 6 Thermal evolution of NSs (A/ = 1.3 Mq) in SXTs from an initial state to their final steady states. Left: evolution, 
at a given mass accretion rate, for several initial values of the inner stellar temperature. Right: evolution from a given initial 
state for several values of the mean mass accretion rate. 



upper limits are clearly much below the predictions of 
the minimal cooling model. In present paper we have 
adopted the same upper l imit on bolometric lum inosity 
of SAX 1808.4-3658 as in lYakovlev et aT 

Thus, despite many theoretical and observational un- 
certainties, it seems that the "minimal cooling model" is 
ruled out by the data on SXTs. Therefore, the mecha- 
nism of the enhanced NS cooling should be more power- 
ful than the process of Cooper pairing of neutrons. 

Nevertheless, mild neutron superfluidity localized in 
the inner NS core is not prohibited. If it is available, to- 
gether with proton superfluidity, the enhanced cooling 
should be of the direct Urea type. If, however, proton 
pairing does not extend to the inner core, neutron su- 
perfluidity remains the only regul ator of cooling o f mas- 
sive NSs. This case was studied bv lGusakov et al" 1(11)05) 
in NS models with open direct Urea process. Because 
of the lack of space, we do not illustrate this case. We 
just note that this model easily fits hot and cold sources 
but scarcely explains quite a large span of intermediate 
sources, for example. Vela and Cen X-4. 



4 Dichotomy of thermal evolution of SXTs? 

As mild neutron superfluidity might exist in the inner 
NS core, let us outline its effect on thermal evolution 
of SXTs. In contrast to INSs, thermal states of NSs in 
SXTs depend not only on NS properties but also on the 
accretion rate. 

As seen from Fig.[5](right), heating curves for medium- 
mass SXTs drop abruptly when the mass accretion rate 
becomes lower than some threshold value specific for an 



NS of a given mass. As we already know, this drop (the 
vertical segment of a curve) signals the onset of neutron 
superfluidity in the NS core when the deep crustal heat- 
ing becomes insufficient to keep the core temperature 
above the maximal critical temperature of neutrons. 

Let us remind that NSs are thermally inertial objects. 
A thermal state of the star as a wh ole may change no - 
ticeably on time scales of ~ 10''' yr IColpi et al.l (|200lh . 
We also recall that heating curves refer to steady-states 
which NSs reach in a few millions years after the on- 
set of the accretion stage in SXTs. The abrupt drops of 
heating curves indicate that thermal evolution of SXTs 
exhibits dichotomy. This means that two NSs with the 
same mass may evolve to very different steady states if 
their mass accretion rates are slightly different. The right 
panel of Fig. [6] illustrates how a change of M by one per- 
cent from some "threshold" value entails a change of the 
steady state thermal luminosity by a factor of ^ 1.5. 

Moreover, a tiny mismatch of masses of two NSs, or 
a small difference of their inner temperatures before the 
onset of the accretion stage, may also induce similar di- 
chotomy. The latter effect is illustrated in the left panel 
of Fig. [6l The former one can be seen in Fig. [5] (and 
Fig. HI): at some "threshold" M, NSs with masses dif- 
ferent by one percent (~ 0.02 Mq) may show thermal 
luminosities which differ by factor of ~ 3. One can also 
speculate that the dichotomy can be caused by a vari- 
able crust composition, which can change the efficiency 
of deep crustal heating at a given M (see Haensel & 
Zdunik, 1990; 2003). Concluding, whatever causes the 
dichotomy, the presence of mild neutron superfiuidity in 
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the NS core may considerably complicate an interpreta- 
tion of thermal emission from SXTs. 



5 Discussion 

We have shown that the data on quiescent thermal emis- 
sion from SXTs require the presence of strong proton su- 
perfluidity and the absence of mild neutron superfluidity 
in the outer NS core (i.e., in low-mass NSs). The mecha- 
nism which controls fast neutrino cooling of massive NSs 
should be more powerful than that due to Cooper pair- 
ing of neutrons and, probably, than the direct-Urca-type 
processes with kaons and pions. We have noticed that 
the presence of mild neutron superfluidity in the inner 
NS core may cause the dichotomy of thermal evolution 
of SXTs. 

Further studies are required to confirm or reject our 
inferences. The data on SXTs are very uncertain. Much 
work is needed to constrain these data and clarify the na- 
ture of variability of some SXTs in quiescence. There are 
still many challenges to respond to in order to build a re- 
alistic theory of deep crustal heating of SXTs. Note that 
current models of accreted crust and deep crustal heat- 
ing are based on one particular cold liquid drop model 
of atomic nuclei and one model of pycnonuclear burn- 
ing. The actual properties of highly neutron-rich nuclei 
surrounded by free neutrons in the inner NS crust, and 
the rate of pycnonuclear fusion of such nuclei (Coulomb 
barrier penetrability, astrophysical 5'-fac tors) are known 
with considerable uncertainty (see, e.g.. lYakovlev et all 
l2006h . In this study we have neglected neutrino emission 
due to Cooper pairing of neutrons in the NS crust in 
order to show directly the effect of proton and neutron 
superfluidity in the NSs core on thermal states of SXTs. 
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